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1 Introduction 

Recall that a Belyi function is a rational function on an algebraic curve with at most 
3 critical values. The corresponding covering of by the algebraic curve ramifies only 
above (at most) 3 points. By fractional-linear transformations, the ramification locus can 
be chosen to be the set {0, 1, oo} c F^. 

According to Belyi [1] and Grothendieck [12], there are deep relations between Belyi 
functions and algebraic curves defined over Q, and dessins d'enfant. 

More generally, one can consider the set of (isomorphism classes of) all coverings on 

with prescribed number of ramified points and with prescribed ramification orders 
above them. Such a topological configuration space is called a Hurwitz space. If we fix 
a ramification pattern for Belyi functions, we typically have a finite set of (isomorphism 
classes of) Belyi functions with the prescribed ramification pattern. If we fix the hypermap 
|32j of 3 permutations for the monodromy group of the covering, the Belyi map is unique. 

In this article we consider coverings of P^ which ramify only above 4 points, and such 
that there is only one simple ramified point in one of the 4 fibers. We refer to these 
coverings as almost Belyi coverings. As is known, Hurwitz spaces for coverings ramified 
only above 4 general points have dimension one \32\ Proposition 3.1]. In fact, any algebraic 
curve can be obtained as some one-dimensional Hurwitz space (with specified monodromy 
permutations) [9j. 

For an almost Belyi covering of degree n, let us denote its ramification pattern by 
-R4(-Pi| -P2I -P3)) where Pi,P2,P3 are 3 partitions of n specifying the ramification orders 
above three points. The fourth partition is assumed to be 2 -|- 1 -1- 1 -|- ... -I- 1. The similar 
notation for a ramification pattern for Belyi maps is i?3(Pi| P2I ^3)) as in [1], |19j . 

The main goal of this paper is to compute generic almost Belyi coverings P^ — > P^ 
with the following ramification patterns: 



i24(5-F5-F5-F5| 2-F2-F2-h2-F2-F2-F2-F2-F2-F2 | 3+3+3-h3-h3-h2-M-hl + l) . (1.3) 

Their degree is 12, 11 and 20, respectively. 

We consider the three specific coverings because of their application to the theory of 
algebraic Painleve VI functions. With certain almost Belyi coverings, one can pull-back 
a hypergeometric differential equation to a parametric isomonodromic Fuchsian equation 
with 4 regular singular points plus one apparent singularity. Equivalently, one may obtain 
isomonodromic 2x2 matrix Fuchsian systems with 4 regular singular points. The corre- 
sponding Painleve VI solutions are algebraic. Knowing suitable almost Belyi maps, one 
can construct explicit examples of algebraic Painleve VI solutions [18] , [2] , [TO] , and solve 
the corresponding isomonodromic Fuchsian equations explicitly in terms of hypergeomet- 
ric functions. More generally, explicit knowledge of any Hurwitz space can be similarly 
used to solve explicitly many types of Fuchsian systems, such as Garnier systems |17| . 

With the three almost Belyi coverings we construct, we pull-back hypergeometric equa- 
tions with the icosahedral monodromy group to isomonodromic 2x2 Fuchsian systems 
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with 4 regular singular points, and the same monodromy group. In total, we compute 
five corresponding algebraic Painleve VI solutions. They have type 37, 38, 41, 42, 43 
in Boalch's classification of icosahedral Painleve VI equations. Three algebraic Painleve 
VI solutions (of type 38, 42, 41, respectively) can be constructed immediatel}0 from the 
almost Belyi coverings with the ramification patterns (jl.ip - ()1.3p . To obtain other two 
algebraic Painleve VI solutions, we compute properly pull-backed corresponding Fuchsian 
systems explicitly. The type 41 solution is related to the Great Dodecahedron Solution of 
Dubrovin-Mazzocco [TT] via an Okamoto transformation. 

Efficient computation of highly ramified coverings or Hurwitz spaces are important 
problems in other fields as well. Therefore these problems attract attention of researchers. 
In [8] , a method is presented to compute one-dimensional Hurwitz spaces (for almost Belyi 
maps, for example) based on the degenerations when 4 ramification loci coalesce into 3 
ramified points. In [22], a computer algebra package is presented for computing genera 
and monodromy groups of Hurwitz spaces or coverings. 

In the next section, we present our computational method. It is basically the same 
method as described in |29l Section 3]. Compared with the most straightforward method 
with undetermined coefficients, we derive equations of smaller degree in undetermined co- 
efficients by using properties of the derivatives of Belyi or almost Belyi maps. We found out 
that very much the same computational method was used in [15], for deriving several rather 
simple Belyi coverings by hand. In Sections [3] through [5] we present our computations of 
the coverings with ramification patterns (|l.ip - (jl.3p . Section [6] demonstrates application 
of the three coverings to computation of algebraic Painleve VI solutions. The coverings 
are i?-parts of -RS-pullback transformations of hypergeometric differential equations to 
isomonodromic 2x2 Fuchsian systems. Computation of i?5-pullback transformations of 
isomonodromic systems 2x2 Fuchsian systems is discussed thoroughly in |27j . 

The authors prepared a Maple 9.0 worksheet supplementing this article, with the formu- 
las in Maple input format, and demonstration of key computations. Readers may contact 
the authors, or search a current website of the first author on the internet, to access the 
worksheet. 

2 The computational method 

Here we briefly recall the straightforward method for computation of almost Belyi coverings 
from to P^, and present an improved method that uses differentiation. To distinguish 
the two projective curves, we write the coverings as P^ — ^ where x and z denote the 
rational parameters of the projective lines above and below, respectively. We assume that 
the three ramification loci indicated in the i?4-notation are z = 0, z = 1, z = oo, in this 
order. We refer to the simple ramification point in the fourth ramified fiber as the extra 
ramification point. 

^Our original motivation for this work was to compute a few missing examples in early versions of [5] of 
icosahedral Painleve VI functions. We did our computations for type 38, 41 solutions unaware of the sixth 
electronic version of [5]. Before the next version of with type 42, 43 examples appeared, we had the 
degree 11 covering and the corresponding type 42 solution as well. Complimentary to [5], computations of 
type 44-45 and 47-52 examples were done independently in [7] and [26j . 
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Let n denote the degree of the covering. By the Hurwitz genus formula [14^ Corollary 
IV. 2. 4], the number of distinct points above {0, l,oo} C must be n + 3 for an almost 
Belyi covering (and n + 2 for a Belyi map); see |18l Proposition 2.1] or [281 Lemma 2.5]. 

The straightforward method to compute an almost Belyi covering with a given rami- 
fication pattern is to write an ansatz of the form 



where F, G, H are general polynomials in x of the factorized form determined by the 
respective partition of n. Specifically, the polynomials have the form Cq YYj=i Pj , where 
C is a constant (undetermined yet), and each Pj is a monic general polynomial of degree 
equal to the number of parts j in the respective partition. Of course, polynomials of degree 
zero can be effectively skipped. To avoid redundancy, one may assume that H is a monic 
polynomial, and may pick 3 of the x-pointj^ as x = oo, x = and x = 1. Expression (j2.ip 
leads to the polynomial identity F = G + H;hy expanding the polynomials and comparing 
the terms to the powers of x one gets a set of polynomial equations. 

This straightforward method was extensively used in [18], [19] to compute Belyi maps 
and almost Belyi coverings of degree up to 12. Those coverings were applied to com- 
pute algebraic transformations of Gauss hypergeometric functions, or compute algebraic 
Painleve VI functions. However, the amount of computations with the straightforward 
method grows quickly for larger n. In particular, the number of variables and algebraic 
degree of initial equations grow linearly with n. The polynomial system may have many 
degenerate (or parasitic ^Q\) solutions, when the rational expression in (12. ip can be sim- 
plified to a rational function of lower degree. Computation of higher degree coverings with 
the straightforward method is hardly possible even with modern computers. In particu- 
lar, our three coverings, including the degree 11 covering, were too hard to compute in 
reasonable time with available PC's. 

Equations of smaller degree for the undetermined coefficients can be obtained by con- 
sidering derivatives of (p{x). According to Couveignes [8], Pricke was probably the first to 
use differentiation to investigate highly ramified maps. More recently, differentiation was 
used for investigation of dessins d'enfant in [2l], [25], [31] and other works. Specifically 
for computational purposes, differentiation was used in [15] and [29] in similar ways. 

A systematic procedure for computation of Belyi maps with differentiation is formu- 
lated in [291 Section 3]. The main trick is to consider logarithmic derivatives of <p and 
ip — 1 from (j2.ip . For instance, the denominator of (p' /(fi is the product of all factors of F 
and H, to the power 1. The numerator is the product of all factors of G, with the powers 
diminished by 1. This gives equations of smaller degree, and easy possibilities for elimi- 
nation. Typically, the degree of equations is diminished by the number of distinct points 
in a corresponding ramified fiber. Compared with computation of Belyi maps, the only 

^Strictly speaking, the x-points for almost Belyi coverings are curves, or one-dimensional branches of 
a generic family, parametrized by an isomonodromy parameter t or other parameter, since the Hurwitz 
spaces for almost Belyi maps are one-dimensional. For simplicity, we ignore the dimensions introduced by 
such parameters, and consider a one-dimensional Hurwitz space as a generic point. 
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adaptation for almost Belyi coverings is that numerators of the logarithmic derivatives 
have an additional degree 1 factor coming from the extra ramified point. 

In this paper, we apply the method in \29\ Section 3] for computation of the almost 
Belyi coverings with the ramification type p.ip - (|1.3|) . In the following Section, we present 
the computational steps specifically for ramification pattern (jl.ip quite in detail. Having 
demonstrated that example, we present our computations for the other two coverings 
with less notice of the routine steps, but we concentrate rather on additional heuristic 
means of solving the obtained systems of equations. In particular, in Section [S] we use 
modular methods for finding coverings with ramification pattern (jl.Sp . Our examples show 
that computational complexity depends not only on the degree of the covering, but also 
on geometric complexity of the solutions (apparently, the geometry is more complicated 
when the degree is prime), or the number of irreducible components of the solutions. 
Accordingly, different heuristic tricks can be useful for different coverings. 

In the rest of this Section, we make a few comments on the fields of definition and 
dimension of Hurwitz spaces. Whether we use the straightforward method or logarithmic 
derivatives, the algebraic equations for the coverings (or Hurwitz spaces) are defined over 
Q. For Belyi functions, the solutions (up to fractional-linear transformations) are isolated 
points, generally defined over an algebraic extension of Q. The field extension may depend 
on the fractional-linear normalization of fixing the points x = oo, x = 0, x = 1. To 
be certain of a minimal Q-extension, one may choose to fix the points of P;J. where the 
ramification order is different from other ramification orders in the same fiber. The central 
question in the theory of dessins d'enfant is how the Galois group of Q/Q acts on Belyi 
coverings (with necessarily the same ramification pattern) or dessins d'enfant. 

General almost Belyi coverings are parameterized by algebraic curves [32], [18] . If 
preferred so, one may consider them as one-dimensional families of almost Belyi coverings. 
The genus of parameterizing curves may depend on the normalization. To get a modelling 
curve of minimal genus, we may strive to fix the points with "isolated" ramification orders 
asa; = oo,a; = 0, x = l, like in the zero-dimensional case. But the equation system can 
be simpler if we adopt the strategy of choosing the points with the highest ramification 
orders as x = oo, x = 0, x = 1. We demonstrate this situation in Section HI 

3 The degree 12 covering 

Here we compute the generic pull-back covering with the ramification type 

i?4(3-h3-h3-hl-M-hl |2-h2-h2-F2-F2-h2|5 + 5 + 2). (3.1) 

Relatively speaking, this is a warm-up example. 

By our conventions, the three partitions of 12 specify the ramification orders above 
z = 0, z = 1 and z = oo, respectively. We choose the simple ramified point above z = oo 
as X = oo, and the simple ramified point above the fourth z-point as x = 0. We do not fix 
X = 1, so there will be the torus action x i— > Ax on the defining equations. The equations 
are expected to be weighted-homogeneous, and the Hurwitz space (if irreducible) should 
be a curve in a weighted-projective space of minimal possible genus. 
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We write the ansatz 



F^G _ „ P 



2 



where 



(/9i2(x) = Co-p-, (pu{x) - l = Goj^, (3.2) 



F = + aix^ + + 03, 



G = + + 62X + 63, 

H = + ClX + C2, 

P = x^ + pix^ + p2a;^ + P3X^ + Pax^ + P5X + P6, 

are polynomials whose roots are the other x-points above z G {0, 1, 00}. In particular, the 
roots of F, H, P are the remaining ramified points. Besides, Gq = liuix^oo ^i2{x) / x'^ is 
an undetermined constant yet. The straightforward method would utilize the following 
consequence of (|3.2p : 

F'^G = P2 + (3.3) 

Following [291 Section 3], we obtain simpler equations in the coefficients of F, G, H, 
P by considering the logarithmic derivatives of ^pi2{x) and ^i2{x) — 1. It is not hard to 
figure out the zeroes and poles of the logarithmic derivatives: 

if 12 FGH (/712 - 1 HP 

Here Ci = C2 = 2 by local considerations at x = cxd. One may generally notice that if 
X = 00 is chosen above z = 00 in a setting like ()3.2p . the constants in the logarithmic 
derivative expressions like in p.4p are equal to the ramification order at x = 00. 
Comparison of the numerators in p.4p gives the following identities: 

2xP = 3F'GH + FG'H - 5FGH', 2xF'^ = 2P'H - 5PH'. (3.5) 



The same type of expressions is derived in computations in [15]. After expanding the 
polynomial expressions and collecting the terms to the powers of x, the first identity gives 
the following equations: 

2pi = 7ci+bi-ai, 

2p2 = 12c2 + 661C1 + 4aici - 2ai6i - 4a2, 

(3.6) 

2pQ = 6016302 + 4026202 - 2026301 + 2036102 - 4a362Ci - 100363, 
= 30263O2 + 036202 - 50363C1. 

The second identity gives the equations 

4oi = 7c 



1, 

4o2 + 2al = \2c2 + 5oipi - 2p2, 



(3.7) 

- 10P6, 
= 202P5 - 50iP6- 



2o| = 402P4 - 30iP5 - 10p6 
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The new equations are sufficient, since tliey are derived from necessary conditions. They 
have smaller algebraic degree, has less degenerate solutions, and can be solved even by 
brute force with Maple's routine solve. More systematically, one may use elimination 
or Grobner basis techniques. The system is overdetermined, but superfluous equations 
only help with Grobner basis computations. As mentioned, the equations are weighted 
homogeneous; specifically 

deg aj = deg bj = deg cj = degpj = j. (3.8) 

The variables pi can be directly eliminated using the first set of equations. This can 
be done similarly for any covering problem with a fiber of only simple ramified points 
(with the ramification order 2) plus possibly one non-ramified point. Also notice that the 
second set of equations does not contain 6j's. There is a dependence between the first two 
equations in (j3.6p and the first two equations in (|3.7p . 

A straightforward way to get the result is the following. Using first equation in (j3.7p . 
we eliminate ai. Then all equations are (still) linear in the ftj's and pi's. We actually get 
12 linearly independent equations in these 9 variables. These variables can be eliminated 
using determinant^ or syzygie^. We get weighted homogeneous equations in ci , C2 , a2 , as . 
We just have to eliminate one more variable! to get an equation for the Hurwitz curve. 
A single such equation in 3 weighted-homogeneous variables is likely to have large degree 
and superfluous factors. To avoid investigating all factors, one may compute two or more 
such equations, and consider only the common factors. It turns out that only the following 
factor gives non-degenerate solutionfl 

IGOalcj + 6912a2C^ - 2256a2C?C2 - 188a2Ci + 103680ci - 81936c? + 20328cfc2 - 1421c?. 

This weighted-homogeneous polynomial defines a curve of genus 0. We can normalize 
ci = 1, and parameterize as follows: 

(2t + 1) (5t + 16) {2t + 5) (15t2 + 2m + 16) ,^ ^, 

= 48^ ' = 16t • ^^-^^ 

Going back, we consequently find parametric expressions for as (and immediately for ai), 
and then for the 9 variables 6i's and pi's. To find the constant Co, we can use (|3.3p 
evaluated at any x £ Q. 



brute way to eliminate the 9 variables is to pick 9 (out of the 12) equations, solve them in the 
9 variables, and substitute into the remaining equations. This is equivalent to computation of 10 x 10 
determinants, with polynomial entries in ci, 02,02,03. In the particular case, this method typically gives 
equations of degree 4 or 5 in 03 alone. 

*Here we mean syzygies between the 12 vectors in the rank 9 free module over Q[ci, C2, 02, 03], with 
the vector components equal to the corresponding coefficients to the 9 variables. When the syzygies are 
applied to the 12 equations, the 9 variables will be eliminated. In the particular case, we can get one 
equation of weighted degree 11 in ci,C2,a2,a3 in this way, and 3 independent equations of degree 12. 

^Here resultants can be used. In our computations, even if we took two equations of degree 4-5 in as, 
Maple 9.5 computed a resultant with respect to 03 in 15-30 seconds. 

^Degenerate solutions are those for which the polynomials F, G, H have multiple or common roots, or 
have the root a; = 0. Consequently, the factors such as 03, 63, C2 or cf — 4c2 can be ignored. It might be 
even useful to search actively for polynomials divisible by the degeneracy factors, so that after dividing 
them out we possibly get polynomials of low degree. In Section[Sl we systematically search for polynomials 
divisible by two resultants defining degeneracy for ramification pattern (|1.3|) . 
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To write the generic solution more compactly, we renormalize x i— > x/4 and multiply 
the polynomials F, G, H by some expressions in t. Here is the covering: 



^u{x) = Cu^^^^, (3.10) 

^12 



where 



~ 16(3t + 4)7 

Fi2 = 2tx^ + Utx^ - 2{2t + + 25t + 16)x - {2t + 5){5t + I6){t'^ + lOt + 6), 

Gi2 = 50t(10t2 + 25t + 16)x^ - 30t{Ut^ - ISi^ - 105t - SO)^^ 

-6(5t + 16)(2t + l)(20t^ + 35*2 + 3t _ iq)x - {2t + lf{f)t + 16)^(5*2 + lOi + 6), 
Hi2 = 3^2 + I2tx + {2t + l)(5i + 16). 



We have 



where 



Pi2 = 20*2(10*2 + 25* + 16)x^ - 12*3(7*2 - 184*^ - 490* - 320)x^ 
+60*2(4068*2 + 4048* + 1024 + 200*^ + 5885*^ + 1386*'^)a;^ 
+20*2(37627*"^ + 100300*2 + 137092*2 + 6415*^ + 250**^ + 92992* + 24576)^2 
-60*(5* + 16) (200*^ + 1260*^ + 3052*^ + 3248*^ + 1126*2 - 137*2 + 256* + 256)2;2 
-30*(5* + 16)2(100*^ + 740*6 + 2289*^ + 3780*"^ + 3600*2 _^ 2040*2 + 700* + 128)x 
-(2* + 1)(5* + 16)2(100*"^ + 740*6 + 2289*^ + 3780*^ + 3600*2 _^ 2040*2 + 700* + 128). 



4 The degree 11 covering 

Here we compute the generic pull-back covering with the ramification type 

7^4(3 + 3 + 3 + 1 + 1 I 2 + 2 + 2 + 2 + 2 + 1 I 5 + 5 + 1). (4.1) 

To get a Hurwitz curve of minimal genus, we can choose the non-ramified point above 
z = 00 as X = 00; the non-ramified point above 2: = 1 as x = 0; and the simple ramified 
point above the fourth z-point as x = 1. Accordingly, we write the ansatz 

ipu{x) = Co^^, ipii{x)-l = Co^, (4.2) 

where F, G, H, P are polynomials of degree 3, 2, 2 and 5, respectively. Zeroes and poles 
of the logarithmic derivatives of ipii{x) and ipii{x) — 1 are easy to figure out, like in ()3.4p . 
The method gives the following identities: 

2{x - 1)P = 3F'GH + FG'H - 5FGH', (4.3) 
2{x - 1)F2 = 2xP'H - bxPH' + PH. (4.4) 
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However, the resulting equations in the coefficients of F, G, H, P are still too complicated 
to solve by direct elimination or Grobner basis techniques. 

Equations of smaller algebraic degree can be obtained if we adopt the strategy to 
normalize the points with highest ramification orders as x = oo, x = 0, etc. (But then 
the Hurwitz space can have non-minimal genus.) Accordingly, we choose the two points 
of ramification order 5 (above z = oo) as x = and x = cxd. We choose the extra ramified 
point above the fourth locus as x = 1, as just above. Then the non-ramified points above 
z = oo and z = 1 are undetermined. We denote their location as x = ci and x = C2, 
respectively. That gives the following ansatz: 

ipu{x) = Go^ -, 99ii(x) - 1 = r, (4.5) 

x^ (x — Clj x^ (x — Clj 

where 

F = x^ + oix^ + a2X + 03, 
G = x2 + 6ix + 62, 

P = + PlX^ + P2X^ + psX^ + P4X + P5, 

and Cq = lima;_^oo ^ii(2;)/2;^ is an undetermined constant. The fractional-linear transfor- 
mation on FI from (I4.5p to (I4.2p is 

Cl(c2 - 1)X + C2(l - Cl) 

X I— > -, — 

(C2 - 1) X + 1 - Cl 

Consideration of logarithmic derivatives of i^ii(x) and (^ii(x) — 1 give the following 
equations: 

5(x - 1)P = 3x(x - ci)F'G + x(x - ci)FG' - (6x - 5ci)FG, (4.7) 
5(x - 1)^2 = 2x(x - ci)(x - C2)P' - (5x2 _ 4^^^ _ g^^^ _^ 5ciC2)P. (4.8) 

Like in the previous example, collecting terms to the powers of x gives a system of equations 
in the Oj's, 6j's, Cj's and pj's. First we pick up the following 8 equations: the terms to 
the powers 5, 4, 2, 1, of x in (14. 7p . and the terms to the powers 6, 1, in (I4.8p . Using 
these equations, we eliminate the a^'s and pj's. Formally, there are 2 solution components, 
but the one with 03 = has to be discarded as degenerat^. We obtain several non- 
homogeneous equations in the 4 variables 61 , 62 , ci j C2 • To get equations for the Hurwitz 

''Here is a stepwise course of elimination. The zeroth powers to x in (|4.7p - (|4.8p . that is, the substitution 
X = 0, gives the equations 5p5 + 5016203 = and Sal = 5ciC2P5. It is easy to eliminate 03,^5. We must 
ignore solutions with 03 = 0, so we are left with aa = — cfc262- Next we consider the coefBcients to the 
first powers to x in (|4.7p - (|4.8|) . and eliminate 02, P4. We get, in particular, 

02 = (5C162 + I2C262 — 3C1C261 — IOC1C2&2) . 

Similarly, we consider the coefficients to the highest degree (5 and 6, respectively) to x in (|4.7|l - H4.8p . and 
eliminate ai,pi. We get, in particular, ai — (36i — 12ci — 5c2 + 10) /II. Having expressed the a^'s and 
Pi, Pi, PS just in terms of &i's and Ci's, the remaining equations are linear in p2, ps. These two pi's can be 
eliminated using discriminants or syzygies, like in the previous Section. 
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curve in two variables, we eliminate bi, 62 by picking up equations of minimal degree in 
them and using resultants. Like in the previous Section, we can compute several resultant 
polynomials (in ci, C2 only) and consider their common divisors as candidate models for the 
Hurwitz curve. It turns out that there is possibly only one component of non-degenerate 
solutions. It is described by a polynomial factor of degree 15 in ci, C2. 

Let Q denote the degree 15 factor. It apparently defines the Hurwitz space for the 
desired almost Belyi map with chosen normalization. The degree of Q in ci alone is 12; the 
degree in C2 is just 6. According to Maple's package algcurves, Q defines a curve of genus 
3. Luckily, the curve is hyperelliptic. Improvised computations produced the following 
Weierstrass model: 

= (St^ + 3t + 2){27f + 71t^ + 130t^ + 140i^ + 120t^ + 64t + 32). (4.9) 

The variables ci and C2 can be parameterized as 

27t^ + 67t^ + 116*4 + 118*3 + 94*2 _^ + 20 + w{3t^ + 2t + 2) 



ci 



C2 



2(t + 2)(t2 + l)(2t2 +3t + 3) 
1107*12 + 7641*11 +Pi _ w{9t^ + 19t'^ + 13t+7){3t^ + 15t^ + 15t^+i5t^+m+26) 



2{t + 2)(3* + l)3(t2 + i)(2*2 + 3t + 3)2(5*2 + 4* + 3) 
where 

Pi = 26055t"+59035tV99475f*+130463t^+138619tVl21015t^+87870f''+51600t^+23798t^+7574t+1460. 
Using these expressions and the equations in bi,b2,ci,C2, we parameterize bi and 62: 
(27*5 -45*^^-190*3 -360*2 -360* -216) ^54^7 _ + w(*+7)(6*2+3*+2)) 



bi 
b2 



32(3* + 1)3(*2 + l)(2*2 + 3* + 3)2(5*2 + 4* + 3) 
P3 + w(* + 7) (6*2 + 3* + 2) (54*"^ - 297*6 - 682*^ - 1 145*^ - 970*^ - 712*2 - 354^ _ ;^q4) 



32(3* + 1)3(*2 + 1)2(2*2 + 3* + 3)(5*2 + 4* + 3)2 
where P2 = 297*^ + 682*^ + 1145*^ + 970*^ + 712*2 _^ 304^ + io4, and 

P3 = 2916t" + 11124t^''^ + 191673t^^ + 764136*" + 1953326t^° + 3445832*^ + 4698345t** + 5040404*'^ 
+4425220*® + 3147872*^ + 1833840*"* + 840864*^ + 301376*^ + 74176* + 11680. 

Here we can stop computations on the hyperelliptic curve. Using fractional-linear trans- 
formation (j4.6p . we can express polynomials H and G in (j4.2p : 

g = x2- (^i-/)(^i+^^) x + ^4^i^, (4.10) 

Cl(C2-l) Cl(C2 - 1)^ 

^ ^ 2 _ (ci-l)(2ciC2+Cibi+C26i+2&2) (Ci-l)2(ci+ 0261+62) ... 

(c2-l)(c2+Ci6i + 62) (c2-l)2(c2+6iCi+62)' ^' ' 

When we write these coefficients in terms of t,w, the square root w conveniently disap- 
pears. Hence the coefficients are just rational functions in *. One may check that the 
algebraic relation between the coefficients of H define an irreducible curve of degree 13 
(and genus 0, as parameterized by *). The projective degree of the parametrization by * 
is 13 as well, so the parametrization is minimal. 
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Now we can get back to the equations induced by (j4.3p - (|4.4p . but knowing parametric 
expressions for 4 variables. Now it is straightforward to find parametric expressions for 
the remaining coefficients of F in ()4.2p . The final expression for ipn can be simplified by 
the renormalization 

+ 15t^ + 15*3 + 45*2 + 40t + 26 

X 7 — 7i n — K ^ X. (4.12) 

(3*2 + 2t + 2)(5t2 +4* + 3) ^ ^ 

Accordingly, the extra ramified point is not fixed to x = 1 in this normalization. The final 
expression can be written as: 

ipu{x) = C^^^^, (4.13) 
-"11 

where 

(2t2 + 3* + 3)(3* + l)2 
108 ' 

77 _ ^3 _ 4941f^+13122t5+19905t'^+17820t3+l0795t2+3962t+879 „2 
-^11 ~ (3t+l)^(2t^+3t+3) 

, (3i+l)(432t^+570t''+330t3-265t2-340f-151) , Q/n+2 , 0+ , Q^^Q+ 1 1 \4 
H 2f-^+3t+3 'JI.^I' -t 01 -\- o) [oT -r I) , 

ri — <'942 I I r)\2^2 I 27t"+270f^+945tg+2160t^+2745t6+1926t5-5f^-1340t3-l440t^-720f-216 ^ 
G-il - [6t +lt + l) X H (2F+3i+3)(3mp ^ 

-4(3* + 1)(2*2 + 3* + 3), 

H^^ = (5*2 + 4* + 3)^2 + 135t«+396tS+715^t^^+790_t3+61^ ^ (2*2 + 3* + 3)(3* + 1)^. 

Let Pii denote the degree 5 polynomial such that 

l-wii(x) = —r^ r^5^. (4.14) 

^^^^ ' 108(2*2 + 3* + 3) Hl^ ^ ' 

We have Pn = (3*2 + 2* + 2) (2*2 + 3* + 3) (3* + l)^^ + . . .. 

As we see, the degree 11 covering is more complicated than the degree 12 covering of 
Section [3l Apparently the geometry of prime degree coverings is more complex. 

5 The degree 20 coverings 

Here we compute generic pull-back coverings with the ramification type 

i?4(5+5+5+5 I 2+2+2+2+2+2+2+2+2+2] 3+3+3+3+3+2 + 1 + 1 + 1). (5.1) 

By fractional-linear transformations, we fix the simple ramified point above z = oo as 
X = oo, and we choose the extra ramified point as a; = 0. 
The ansatz is 

i?5 p2 

(^20(2;) = C'o^g^, ^iQix) - 1 = C'o^g^- (5.2) 

where 

F = rr^ + a\x^ + 02^2 + a^^x + 04, 

G = x^ + + h2X^ + 63x2 + 64X + 65, 

H = x^ ^ c\x^ + C2X + C3, (5-3) 

P = + piX^ + P2X^ + . . . + V'iX + pio, 
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and Co = liniaj^oo '/^(x)/^^ is an undetermined constant. Consideration of logarithmic 
derivatives of ^p2oix) and v?2o(a^) — 1 give the following identities: 

2xP = bF'GH - SFG'H - FGH', 2xF^ = 2P'GH - WG'H - PGH' . (5.4) 

After expanding and collecting terms to the powers of x, we get a system of equations 
in the coefficients of F, G, H, P. Since we do not fix x = 1, the equations are weighted 
homogeneous, with the same grading as formulated in (13. 8|) . 

There are many possibilities to eliminate variables from the equation system. For 
example, one may use the first 10 equations of the first identity to eliminate all pj's. 
Then we can use the first equation from the second group to eliminate bi; the subsequent 
equation turns out to be void. The next 3 equations allow us to eliminate 63, 64 and 65. 
But still, there are too many variables left to solve the system by force. 

Our strategy is the following. We solve the equations modulo several large primes, iso- 
late non-degenerate solutions, and we try to lift them to the characteristic 0. The principle 
aim is to derive modular polynomial equations which characterize only non-degenerate so- 
lutions. When lifted to Q, those equations are expected to have low degree and rather 
small coefficients, because they would contain information only about the relevant solu- 
tions. Eventually, it turns out there are a few connected components of non-degenerate 
solutions; we are able to separate them on the modular level, so that each lifted equa- 
tion system describes only one connected component. We effectively avoid intermediate 
computations with huge Q-coefficients, and consider over Q only those equation systems 
which describe isolated components of the generic solution. Actually, we are able to get 
just a few new polynomial equations over Q of low degree, but that is just enough for a 
breakthrough simplification of the original system. 

We use the computer algebra package Singular [13], well suited for ring-theoretic ma- 
nipulation modulo large primes. For a prime number p, let Fp denote the finite field with 
p elements. We did computations modulo these primes: 

p £ {32003, 31991, 31981, 31973, 31963}. (5.5) 

The solutions we found after considering the first 4 primes, while computations modulo 
31963 were done for checking purposes only. 

We use the weighted grading as formulated in (j3.8p . To be able to discard degenerate 
solutions, let Q denote the resultant of F and G with respect to x, and let Z denote the 
resultant of F and H. Their weighted degree is 20 and 12, respectively. A solution of the 
original equation system is degenerate if and only if Q = or Z = 0. 

For each prime number p from (15. 5p . we do computations in two rings: 

Ri = Fp[ai,a2,a3,a4,62,ci,C2,C3,Z], 
i?2 = Fp[ai,a2,a3,a4,62,ci,C2,C3,Z, Q]. 

We assume that 61,63,64,65 and pi's are eliminated from the original system. Let Jq 
denote the graded ideal in R2 generated by the original polynomial equations (after the 
elimination) and by definitions of Z and Q. We wish to find polynomials in Jq divisible 



12 



by Z or Q, and that we could get equations of lower degree for non-degenerate solutions 
by dividing such polynomials by the factor Z or Q. 

Let J\ denote the restriction of Jo onto -Ri . For the beginning, we compute a Grobner 
basis for J\ in R\ with respect to the total degree reverse lexicographic ordering with 

a/i>- cz>- a'i>- C2>- h2>- a2>- c\>- axy- Z. (5.6) 

With this ordering, a (weighted) homogeneous polynomial is divisible by Z if and only if 
the leading term is divisible by Z. Let G\ denote the Grobner basis for J\. Actually, we 
computed G\ up to bounded degree 25. The computations were done on a Dell laptop 
computer with Pentium M 1700MHz processor on the Windows XP platform. 

The first element of G\ divisible by Z occurs in degree 24. (Computations up this 
degree take 430 seconds.) Non-degenerate solutions should satisfy the other degree 12 
factor. We have two options: either use the new degree 12 equation immediately and 
recompute the Grobner basis through degrees 12 to 24, or continue computations in degree 
25 in the hope of finding more elements of G\ divisible by Z . The second option appears to 
be more acceptable since its next step takes less time (510 versus 650 seconds). Besides, 
the more greedy strategy of using lowest degree new polynomials immediately leads to 
more frequent and lengthier recomputations of Grobner bases. In general, one may try 
different tactical choices when making computations modulo first few different primes, and 
then use the best options when computing modulo other primes. 

Our computations are summarized in Table [TJ Recall that the Hilbert series of a graded 
ring R is the series X^^o where hj is the dimension (over the ground field) of the jth 
graded part of R. We refer to the numbers hj as Hilbert dimensions. 

In the first column of Table [H we list the weighted degrees from 11 to 25. The 
numbers in small font are the Hilbert dimensions for R\ (i.e., the number of monomials 
of the weighted degrees in i?i.) In the second column, we give the number of elements of 
G\ up to each degree, and (in the small font) the Hilbert dimensions for the ring R\l J\. 
As mentioned, there is one Z-multiple in degree 24 and two independent Z-multiples on 
degree 25; this is indicated in the next-to-last row, and by the bold face of the numbers 
above. Therefore we have 3 new suitable equations of degree 12 and 13. Let Ji denote 
the ideal generated by J\ and the 3 new equations. 

Subsequently, we compute the Grobner basis for J2 up to degree 24. Column 3 of 
Table [T] gives the same size statistics for J2 . This second run gives 2 new equations of 
degree 12. We iterate the procedure of adjoining new equations and recomputing Grobner 
basis until we don't see Grobner basis elements with the leading monomial divisible by 
Z. As indicated by Table [H the subsequent two runs give bonanzas of 83 and 101 new 
equations of degree 10 to 12. The same quantity and degree of new equations occurs 
modulo each chosen prime, which is a good indication. The statistics of Table [1] indicate 
the complexity of computations in each run. In particular, the Hilbert dimensions indicate 
the size of Grobner basis elements in each degree. 

After the 6th run, we don't get Grobner basis elements divisible by Z . Then we redo 
this run in the ring i?2; we use the same ordering as in (j5.6p with additionally Z y Q. 
With this ordering, a homogeneous polynomial is divisible by Q if and only if the leading 



13 



J— /C^i CC- 


Grobner basis (re)computations 
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317 




17827 
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363 
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1415 
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852 
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3645 
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2 


83 


4+97 


9+3 


1+2 


7+5+1 


Time (s) 


941 


918 


1367 


520 


142 


93 


2 



Table 1: Classical transformations of hyperbolic hypergeometric functions 
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term is divisible by Q. If the degree of the polynomial is less than deg Q = 20, we still 
have the same criterium for divisibility by Z. We do the computation in R2 up to degree 
24, and get one Grobner basis element of degree 23 and two elements of degree 24 divisible 
by Q. After dropping the factor Q, we get an equation of degree 3 and two equations of 
degree 4 for non-degenerate solutions! We feed these equations into the current Grobner 
basis in R2; the following run gives a few equations divisible by Z of degree up to 19. 
After adding the new equations (of degree up to 19 — 12 = 7), we get a stable Grobner 
basis: 17 equations of degree up to 7, plus expressions for Z and Q. This looks like the 
wanted polynomial equations for non-degenerate solutions. 

The computations with Singular were semiautomatic. For the first two primes, we tried 
a couple of strategies when to recompute a Grobner basis. For computations modulo other 
primes we chose the most effective degree bounds. Those computations were automatic 
in principle; in particular, the 6th run was done directly in R2. As mentioned, the new 
equations in corresponding recomputations have the same degree (and leading monomials) 
modulo each chosen prime. Computations for one prime take about an hour. 

We would like to lift the simplest new equations to characteristic 0. Modular lifting of 
Grobner bases is considered by several authors; see for example [3] and further references. 
Strictly speaking, we do not know good bounds for the size of the Q-coefficients. Sure 
bounds must be huge, due to suspected complexity of intermediate computations purely 
in Q. Our whole idea is to escape the intermediate computations by considering over 
Q exclusively equations for non-degenerate solutions only. Eventually, we get actually 
correct solutions after using just the few primes and a few new lifted equations. At the 
end of this Section we indicate a way to check for certain that there are no other solutions. 

A straightforward way to lift equations to characteristic is to lift the coefficients to 
rational numbers with smallest numerators and denominators. Modular reconstruction of 
rational numbers is a well-known problem; the basic algorithm is proved in [30]. A good 
indication that a lift of a polynomial is correct is that the denominators of the coefficients 
have the same factors. If the straightforward method did not work, we tried the LLL 
reduction algorithm [21] for the lattice generated by: the vector whose entries are the 
integer coefficients modulo the composite modulus, plus one extra zero component; and 
the vectors with two non-zero entries — the entry 1 at the extra component, and the 
composite modulus at some (and each) other component. In successful cases, the shortest 
LLL basis vector was shorter by several orders of magnitude than the other vectors, which 
is a convincing indication. We used the extra vector component because Maple's LLL 
routine requires linearly independent input. 

As mentioned, each 6th run of modular computations gives one equation of (weighted) 
degree 3 and two equations of degree 4. The degree 3 polynomials can be lifted quite 
easily. Up to a scalar factor, the result is convincingly 

625 „ 875„ I 11„ ^ I 99A ^ 643 _ „ 99 „3 179 _ „ i 44„ t 

— C3 ^03 + iiC2Ci + II02C1 ^a2Ci — yqC^ ^aiC2 + -3-0102 

_l369 _ _ I 143„ ^2 _|_ 3811 „2„ 895 „3 

The two polynomials of degree 4 can be reliably lifted from the obtained modular data as 
well. With the lifted polynomials, we can immediately eliminate C3 and 04. But further 
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brute force computations seem to be too cumbersome still. 

The final Grobner basis for non-degenerate solutions appears to have new equations of 
degree 5 or higher, but the modular data is not sufficient to lift them. We proceeded then to 
eliminate variables modulo the primes, so to get modular equations in a minimal number 
of variables. These equations indicated that there are a few irreducible components of 
non-degenerate solutions. It was also reasonable to expect that these equations would be 
easier to lift, even if of higher degree. 

Specifically, we used Grobner basis computations with respect to an elimination or- 
dering, still in Singular. The Hurwitz space should be a curve (possibly reducible), so we 
must have weighted homogeneous equations in three variables. The equation in oi, ci, 02 
has degree 14 modulo each prime in (jS.Sp . It factors as follows: two factors of degree 8 
and 6 modulo 31991 (and 31963); or three factors of degree 8, 3, 3 modulo other primes. 
This suggests that there are three families of non-degenerate solutions, two of them are 
conjugate over a quadratic extension of Q. It was possible to lift the degree 8 and degree 
6 = 3 -|- 3 factors to Q using the LLL algorithm. 

The degree 6 polynomial factors over Q(\/— 15) into two factors of degree 3, like ex- 
pected. One may check that the factors define genus curves. They can be parameterized, 
and the other variables can be uniquely parameterized as well using the original equations 
and the three lifted equations of degree 3 or 4. But one may notice a shortcut: ramification 
type (j5.ip can be realized by a composition of two coverings of the types: 

i?3(3-hT + T|5| 2-^2-^1) and R^iS + 1\2 + 1 + 1\2 + 2). (5.8) 

The hats in the first expression indicate the ramification locus of the subsequent degree 4 
covering. The coverings of these two types are known due to their application to Gauss 
hypergeometric functions and algebraic Painleve VI functions, as we recall in Section [6] 
below. The (normalized generic) coverings for these ramification types are: 

(5-3v^) (128x + 7 + 33v^)^ 

Vsl^J = ; , , (5.9) 

8000x(10242;-781-17V-15)'^ 



(t-3)3(3t- l)3(z + l)(z + t) 

(t-iy {{t-iyz + t{t+i)) (4z + 3{t+i)y 



"P^y^l = -xo... -xo — . ... ...... — T^TTT— TTTs- (^-^O) 



Note the appearance of \J —\h in the first covering. The second covering has the parameter 
i as it is an almost Belyi function. The composition 994 o 995 (x) can be written in the form 
(j5.2p . One may check that the coefficients ai,ci,a2 in this form parameterize the degree 
6 polynomial factor (or one of the degree 3 factors). It follows that the two coverings 
implied by the degree 6 factor are (/94 o (^5 and the conjugated version. 

The degree 8 factor in ai, ci, defines a genus curve as well. Its parametrization 
gives rise (by the original equations and the three lifted equations of degree 3 or 4) to the 



16 



following impressive solution of the covering problem: 



ai = 4(3t^ + 9t3 + 43t2 + 40t + 12), (5.11) 

02 = 6(182**^ + 728i^ + 2373*"^ + 3584*^ + 2632^^ + 960t + 144), 

ag = 4(i + l)(8t2 + 7t + 2)(1029t^ + 3246t^ + 9608t^ + 10224*^ + 4752t + 864), 

a4 = {5t + 6)(8t2 + 7t + 2)2(1029t^ + 3246*'' + 9608t^ + 10224*^ + 4752t + 864), 

bi = -i(4i4 + 12*3 _ 859^2 _ 892t - 276), 

62 = -i(448t^ + 1756*5 - 59175*^ - 123710*^ - 101900*^ - 39144* - 6048), 

h = -i(8*2 + 7* + 2) 

x(2352*^ + 9399*5 _ 256240*^ - 577360*^ - 506160*^ - 204336* - 32832), 

64 = -|(8*2 + 7* + 2)2(56*2 + 79* + 34) (49*^ + 132*^ - 5184*^ - 4752* - 1296), 

65 = -|(8*2 + 7* + 2)3(49*2 + 76* + 36) (49*^ + 132*^ - 5184*^ - 4752* - 1296), 
ci = 1(100*^ + 300*3 + 507*2 +388* + 108), 

C2 = i(5* + 6)(8*2 + 7* + 2)(140*3+307*2 +308* + 108), 

C3 = i (5* + 6)2 (8*2 + 7* + 2)2 (49*2 + 76* + 36). 

From now on, let us denote by ip2o{x) the covering defined by (j5.2p . (j5.3p and ()5.1ip . We 
refer to the other two solutions as compositions of the degree 4 and 5 coverings. 

We looked for the degree 20 coverings by modular methods, deliberately ignoring the 
size of intermediate would be computations over Q. It is theoretically possible that some 
generic solutions are missing, since they coincide with the derived solutions modulo each 
of the considered primes, or the original equation system is insufficient modulo those 
primes. This possibility has extremely low probability. Very likely, the "bad" primes can 
be only a few small prime numbers. Regarding the application to Painleve VI functions 
of the next section, we already know that there are exactly two algebraic solutions of 
Pvi{0, 0, 0, —2/3; *) from the work of Dubrovin-Mazzocco [11]: the Cube solution and the 
Great Dodecahedron solution. After an Okamoto transformation we have exactly two 
solutions of Pi//(l/3, 1/3, 1/3, 1/3; *); they are obtainable from the composition 934 o ip5{x) 
and the irreducible covering ip2o{x), as we will see. 

To check for certain that there are no other covering with the ramification pattern (j5.ip , 
the canonical method is combinatorial. Each branch of the Hurwitz space corresponds to a 
4-tuple of permuations of 20 elements, of cycle types 5+5+5+5, 3+3+3+3+3+2+1+1+1, 
etc. The braid group on 4 braids acts on the branches of the same connected component. 
There must be only three orbits of the braid group, giving three connected components of 
the Hurwitz space, corresponding to the composition 994 o 995(2;), its complex conjugate, 
and Lp2o{x). This method is strict [8], [22], but it requires computation of all permutation 
combinations with the given cycle type and identity product. More geometrically, one 
may introduce deformations of dessins d 'enfant [18] , [19] , cacti [23] or similar geometric 
objects [6} pg. 105] that represent almost Belyi coverings in the same way as usual dessins 
d'enfant correspond to Belyi maps, observe homotopic action of the braid group, and count 
possible "deformation" drawings with the given branching type. 
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6 Application to algebraic Painleve VI functions 



As noticed in [28], [29] and [18], [T], certain Belyi coverings occur with algebraic trans- 
formations of Gauss hyper geometric solutions. These transformations are induced by 
pull-back transformation of a hyper geometric differential equation to a hypergeometric 
equation again. In particular, Belyi covering (j5.9p transforms between standard hyperge- 
ometric equations with the icosahedral and tetrahedral monodromy groups. Here is an 
induced hypergeometric identity: 



2F1 



1/4,-1/12 _ ^ 7-33^^^^ ^'^'^ ^ 11/60, -1/60 
2/3 



IH x\ 2F1I 

128 ) V 2/3 



1 



(6.1) 



This is the same transformation as formula (50) in [28], but with a different definition of 
935(2;). This formula can be checked by comparing the Taylor expansions of both sides 
around x = 0. 

Similarly [T5], [TU], almost Belyi coverings with certain ramification patterns can be 
used to pullback hypergeometric differential equations to 2 x 2 isomonodromic Fuchsian 
systems with four singularities. Correspondingly, one may derive algebraic solutions y{T) 
of the sixth Painleve equation: 

(fy 1/1 1 I \ ( dyV f I 1 I \ dy 



dT^ 2\y^y-l^y-T)\dT) [t ^ T - I ^ y - T J dT 

y(y-l)(y-T) / T T-l T(T - l)\ , , 

where a, /9, 7, (5 G C are parameters. The standard correspondence between solutions of 
the sixth Painleve equation and the mentioned isomonodromic Fuchsian systems is due to 
Jimbo-Miwa [16]. If the singular points of the Fuchsian system are x = 0, x = 1, x = T, 
X = 00, and the local monodromy differences at them are, respectively, ^ij ^t-, Ooo, then 
the corresponding Painleve equation has the parameters 

« = ^ , /3 = -y, 7 = y, ^ = (6-3) 

We denote the corresponding Painleve VI equation by Pvi{Oo,9i,6t,0oo',T). 

General pullback transformations of 2 x 2 Fuchsian systems d^{z)/dz = M{z)'^{z) 
have the following form: 

z^R{x), -^{z)^ S{x)-^{R{x)), (6.4) 

where R{x) is a rational function of x, and ^(x) is a Schlesinger transformation, usually de- 
signed to remove apparent singularities. For transformations to parametric isomonodromic 
equations, R{x) and ^(x) may depend algebraically on parameter(s) as well. In [T7], [15] . 
|27| . these pullback transformations are called RS -pullback transformations, meaning that 
they are compositions of a rational change of the independent variable z 1— > i?(x) and the 
Schlesinger transformation S{x). The Schlesinger transformation ^(x) is analogous here to 
projective equivalence transformations y(x) 9{x)y{x) of ordinary differential equations. 
If ^(x) is the identity transformation, we have a direct pullback of a Fuchsian equation. 
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If z = R{x) is an almost Belyi covering with a suitable ramification pattern, one can 
pick up hyper geometric equations (in a correspondingly normalized matrix 2x2 form) 
and choose appropriate Schlesinger transformations S{x) so that the pullbacked Fuchsian 
equation would be isomonodromic and have four singular points, and there would be a 
corresponding algebraic solution of the sixth Painleve equation. These iJS-transformations 
are defined in [18], [19]; their algorithmic construction is considered thoroughly in |27] . 
The notation for suitable classes of these i?S'-pullback transformations is 



f eo 


ei 


Coo 




Pi 


P 

^ oo 



(6.5) 



Here the subscripts 2 and 4 indicate a second order Fuchsian system with 4 singular 
points after the -RS-pullback; Pq, Pi, P^o define the ramification pattern R4,{Po \ Pi \ Poo) 
of the almost Belyi covering R{x); and eo, ei, eoo are the local exponent difference of the 
hypergeometric equation. 

With the almost Belyi coverings 9912 (x), ipii{x), (p2oix) of this paper, we can construct 
-RS-transformations of the types 



n 02 / 1/3 
002/ 1/3 

Jt04 I 3_|_3_,_3_^;^_,_;^ 



1/2 

2+2+2+2+2+2 
1/2 

2+2+2+2+2+1 



1/5 
5+5+2 / ' 



1/5 
5+5+1 



r.c-2 / 1/5 
it>J4 I 5_,_5_,_5_,_5 



OC.2 / 1/3 

3+3+3+1+1+1 

r>a2 ( 1/3 

^^3_|_3_^3_|_;^_,_;^ 

1/2 



1/2 

2+2+2+2+2+2 
1/2 

2+2+2+2+2+1 



2/5 
5+5+2 

2/5 
5+5+1 



2+2+2+2+2+2+2+2+2+2 

and derive algebraic solutions of, respectively. 



1/3 

3+3+3+3+3+2+1+1+1 ; ' 



Py/(l/3,l/3,l/3,3/5;r), 
Py/(l/3,l/3,l/2,4/5;T), 



Py/(1/3,1/3,1/3,1/5;T), 
Pv/(l/3,l/3,l/2,2/5;r), 
Py/(1/3,1/3,1/3,1/3;T). 



(6.6) 
(6.7) 
(6.8) 



All of the i?5-pullbacks transform hypergeometric equations with the icosahedral mon- 
odromy group to isomonodromic Fuchsian systems with four singular points and the same 
monodromy group. The Painleve VI solutions are called icosahedral [5]; there are 52 types 
of them up to branching representation of the icosahedral monodromy group, or Okamoto 
transformations. The solutions of (I6.6p - (l6.8p have following Boalch types, respectively: 38, 
37, 42, 43, 41. As mentioned in a footnote to the introduction section, our computations 
of these solutions by the method of iJS-transformations is independent from [5] . 

Direct results relating i?5-pullback transformations to algebraic Painleve VI solutions 
are presented in [27]. The most conveninet results are reproduced here. 

Theorem 6.1 Let kQ,ki,koo denote three integers, all > 2. Let : F^. ^ denote an 
almost Belyi map, dependent on a parameter T . Suppose that the following conditions are 
satisfied: 

(i) The covering z = (p{x) is ramified above the points z = 0, z = 1, z = 00; there is 
one simply ramified point x = y above P?, \{0, 1, 00}; and there are no other ramified 
points. 
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(a) The points x = 0, x = 1, x = oo, x = T lie above the set {0, 1, oo} C P^. 

(Hi) The points in ip~^{0) \ {0, l,r, oo} are all ramified with the order /cq. The points 
in ip~^{l) \ {0, 1,T, oo} are all ramified with the order ki. The points in ip~^(oo) \ 
{0, 1,T, oo} are all ramified with the order k^o- 

Let ao,ai,aT)floo denote the ramification orders at x = 0, 1,T, oo, respectively. Then the 
point X = y, as a function of x = T, is an algebraic solution of 

P f °0 Ql QT -i QoO rj.\ fr, o^ 

^yi\T. — 'I — 'I — — 5^ • ^^-^^ 

Proof. This is Theorem 3.1 in |27j . □ 

Theorem 6.2 Let z = ip{x) denote a rational covering, and let F[x), G{x), H[x) denote 
polynomials in x. Let E denote the hypergeometric equation with the local exponent dif- 
ferences Co, ei, Coo at, respectively, z = 0, z = I, z = oo, Suppose that the direct pullback 
of E with respect to ^{x) is a Fuchsian equation with the following singularities: 

• Four singularities are x = 0, x = l,x = oo and x = T, with the local monodromy 
differences do, di, dx, doo, respectively. The point x = oo lies above z = oo. 

• All other singularities in \ {0, 1,T, oo} are apparent singularities. The apparent 
singularities above z = (respectively, above z = 1, z = oo) are the roots of F{x) = 
(respectively, of G{x) = 0, H{x) = Oj. Their local monodromy differences are equal 
to the multiplicities of those roots. 

Let 5 denote a non-negative integer such that A := deg-F + degG + degif + 5 is even. 
Suppose that {U2,V2,W2) is a syzygy between the three polynomials F, G, H, satisfying, 
if 6 = 0, 

degC/2 = 4 -degF, deg I/2 = # - deg G, deg W2 < f - deg if, (6.10) 
or, if 6 > 0, 

deg C/2 < ^ - deg F, deg V2 < ^ - deg G, deg W2 = ^ - deg H. (6.11) 
Then the numerator of the (simplified) rational function 

U2W2 f{eo - ei + eoo) {FU2)' ^ {HW2)' \ ^ (cq - ei - e^,) ^21^2 ^' 



G \ 2 if FU2 HW2 J 2 F if-l 

(ep + ei - Coo) U2V2 ^' 

+ 2 ir^(^y 

has degree 1 in x, and the x-root of it is an algebraic solution of Pyiido, di, dx, doo -\-6;t). 

Proof. This is Theorem 6.1 in |27| . □ 

The first theorem here is actuahy a special case of the second one, when the syzygy 
(U2, V2, W2) has one of the components equal to zero. The Painleve VI solutions obtained 
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l/fcl 1/feo 
Pi Poo 

solutions can 



5+5+2 ) '^ith respect 



in this case can be seen as inverse functions of particular projections [HI Figure 5] of the re- 
spective Hurwitz spaces to f\ . The implied i?S'-transformation is RSl ^ ^^^^ 
where Pq, Pi, Poo define the ramification pattern of (p{x). The Painleve V 
be derived without computing the Schlesinger part S{x) of the i?5-transformation. It is 
these solutions that are implied or computed in and [18] . 

The more general Theorem 16.21 means that the same almost Belyi covering can be 
used to pullback several hypergeometric equations, and hence derive several algebraic 
Painleve VI equations. For each implied i2S'-transformation, the Schlesinger transforma- 
tions depend actually on two syzygies for the same polynomial triple {F, G, H), but a single 
Painleve solution depends on one syzygy, as stated in the theorem. (The other syzygy can 
be usually used to compute an algebraic solution of other Painleve VI equation.) 

Both theorems have to be applied to fractional-linear normalizations of the coverings 
^uix), 9911 (x), ip2oix), where three (of the four) singular points of the transformed Fuch- 
sian equation are chosen to be rc = 0, a; = 1, x = oo. Theorem 16 . 1 1 eventually gives solutions 
of Py7(l/3,l/3,l/3,3/5;r), Py/(l/3, 1/3, 1/2, 4/5; T), Py^l/S, 1/3, 1/3, 1/3; T). These 
are icosahedral solutions of Boalch types 38, 42, 41, respectively. 

In particular, an i?S'-pullback RS^ (^^^^^^/^i^i^i 2+2+2+2+2 
to 2; = 9312(2;) is a Fuchsian system with singularities at a; = 00 and the roots of Gi2{x) = 0. 
The local monodromy differences are, respectively, 3/5, 1/3, 1/3, 1/3. A suitable normal- 
izing fractional- linear transformation may leave x = 00 invariant, but it must move two 
roots of ^12(2;) to the locations x = and x = 1. Finding one root of G12 is equivalent to 
considering Gi2{x,t) = as an equation for an algebraic curve. The curve has genus 0; it 
can be parametrized as follows: 

^ _ 3{2s + 1) (8g^ -2s + 17)(2g^ + 2g + 3)(g + l)^ 

2(s3 + 3s + l)' ^~ 10(s3 +3s + l)(2s2 + S + 2) ' ^' ' 

After the repameterization of t by s, the other two roots of Gi2{x) are equal to 

_ (s-2)(4s+l)(32sfi -8s5+164s-'-94g3+9ig2^2s+18) i {s-2)^ (is+l)^ {8s^ +6s-l) w 
^ ~ 20{s^+3s+l)(4s^-s+l)^ 20{2s+l)(s^+3s+l){4s^-s+l)^ ' 

where w = a/(s - 2)(2s + l)(2s2 + s + 2). It appears that the polynomial ^12(2;) splits 
over the function field of a genus 1 curve. Let us denote the x-root in (j6.13p by cq, and 
the latter 2 roots by c+, c_. A normalizing projective coordinate for is: 

Ai2(x) = ^^fiL. (6.14) 

c_ — c+ 

Equivalently, the normalizing fractional-linear substitution to the new coordinate is given 
by Aj~2^(x): x c+(l — x) + c_x. By Theorem 16.11 applied to 99i2(Aj~2^(x)), an algebraic 
solution y38(T38) of Py/(l/3, 1/3, l/3,3/5;r38) is given by 

^38 = Ai2(co), yss = Ai2(0). (6.15) 
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Explicitly, we have 

1 3(2s + l)(32s'^ + 32s6 + 138s5 + 25s'' + 130s3 + 30s2 + 20s-10) 
T'is, — — -| , — , 

2 2(s - 2)2(4s + l)2(2s2 + s + 2) V(s - 2)(2s + l){2s^ + s + 2) 

1 (2s + l)(32s6 - 8s^ + 164s^ - 945^ + Ols^ + 2s + 18) 

2/38 = t: + 



2 2(s - 2)(4s + l)(8s3 + 6s - 1)^/(5 - 2)(2s + l)(2s2 + s + 2) 

To get the parametrization obtained in [5], one has to change s — > — (s + l)/2s and the 
branch of the square root. 

Similarly, to apply Theorem 16. II to the degree 11 covering ipii{x), we have to compose 
it with a fractional-linear transformation of F^. which leaves x = oo invariant, and moves 
the roots of Gii{x) from (j3.10p to the locations x = and x = 1. The roots of Gii{x) are: 

27*1° + 270*9 + 945*^ + 2160t^ + 2745*^ + 1926*^ - bt^ - ISAOt^ - 1440*^ - 720t - 2 16 
^ ~ 2(3t + 1)2(3*2 + 2t + 2)2(2*2 + 3t + 3) 



_^3(t + 2)2(t2 + 1)2(3^2 + 2t + 2) Y^3(t + 2)(t + 7) (3*2 + 3t + 2) 
2(3t + 1)2(3*2 + 2t + 2)2(2*2 + 3t + 3) ' 

Let us denote these two roots by c+ and c_. Then we can use the same expression 
(j6.14p for the normalizing fractional-linear transformation. To distinguish, we denote this 
fractional-linear transformation by Aii(x). Now, let yo denote the extra ramification point 
of 9311, outside the fiber of {0, 1, oo} C P^. We have: 

_ 3t^ + 15t^ + 15*3 + 45*2 + 40* + 26 
^° ~ (3*2 + 2* + 2)(5*2 + 4* + 3) ■ 

By Theorem l6.1l applied to (^ii(Aj~/(x)), an algebraic solution of Pi//(l/3, 1/3, 1/2, 4/5; T42) 
is given by 

T42 = Aii(0), y42 = Aii(yo), 
We arrive at the following Painleve VI solution: 

^ _ 1 27*i° + 270*9+945*'^ + 2160*^+2745*6-hl926*5-5*^-1340*3-1440*2-720*-216 

~ 2 ^ 6(* + 2)2(*2 + l)2(3t2 + 3i + 2) V3(* + 2)(* + 7)(3*2 +3tT^ ' 

1 (* + 7) (45*6 + 144*5 _^ 258*^ + 228*^ + 121*2 _^ 24* - 12) 
y42 = - + 



2 6(5*2 + 4t + 3)(t2 + l)(t + 2) Y^3(* + 2)(* + 7) (3*2 +3t + 2) 

The solution has genus 1 as well. To get the parametrization of the same solution obtained 
in [5], one may substitute * —(2s — l)/(s + 2). 

To get a solution of Py/(l/3, 1/3, 1/3, 1/3; T) by Theorem 16. 11 the covering (p2o{x) can 
be composed with a fractional-linear transformation of Fl. which leaves x = 00 invariant, 
and moves two roots of H2o{x) to the locations x = and x = 1. Like in the case with 
99i2(x), one root of H2o{x) can be made explicit by parametrizing the curve H2o{x,t) = 0: 

^_ 2(2s3+4s2-4s + 3) 16(s-2)2(2s2 + s + 2)2(s2-2s + 6) 

5(2s-l)2 ' 75(2s-l)4 • ^^-^^^ 
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After the repameterization of t by s, the other two roots of H2q{x) are equal to 

16s(8g^ - lis + 8)(56g'^ - ims'^ + 318g^ - 269s^ + 31g^ + 75^^ - 285 
^ ~ 25(2s - 1)8 



144s(s - l)(8s2 - lis + 8)(2s3 + 45^ - 4s + 3)^s(8s2 - lis + 8) 
^ 25(2s- 1)8 ' 

Let us denote the x-root in (I6.16P by cq , and the latter 2 roots by c+ , c_ . Then a suitable 
projective parameter \2{){x) is given by the same expression as on the right-hand side of 
(j6.14p . An algebraic solution y4i(T4i) of Pyi{l/?,, 1/3, 1/3, 1/3; T41) is given by 

?4i = A2o(co), 2/41 = A2o(0). (6.17) 

Explicitly, we have 

1 (s + l)(32s8-320s'^ + 1112s6-2420s5+3167s^-2420s3 + 1112s2-320s+32) 
^41 — t: "I : ~ — , . „ . . 3 ) 



2 54s(s - l)(v^s(8s2 - lis + 8))' 

1 _ 56s^ - 166s'^ + 318s^ - 2695^^ + 31s^ + 755^ - 28s + 8 

2 18s(s - l)(2s3 + 4s2 - 4s + 3) Vs(8s2 - lis + 8) 



To get the parametrization presented in [5l Theorem C], one has to change s — > 1/s and 
the branch of the square root. This solution is related via an Okamoto transformation 
to Great Dodecahedron Solution |llt pages 134-143]. The Dubrovin-Mazzocco solution 
solves Pv'/(0,0,0,-2/3;r4i). 

The two composite coverings (^4 o (^5 of degree 20, with the same ramification pattern 
as ip2o, generate the same algebraic Painleve VI solution as the degree 4 covering 994 de- 
fined in (jS.lOp . That algebraic solution of the same equation Py7(l/3, 1/3, 1/3, 1/3; T41) 
is parametrized in [181 Section 3.2], basically using Theorem 16.11 Via the same Okamoto 
transformation, we get the Cube solution in [11], of Py/(0, 0,0, — 2/3;T4i) as well. Com- 
position of an almost Belyi covering with a Belyi covering never changes the algebraic 
Painleve VI solution. Complimenting computations in [181 Sections 3.1, 3.4], up to 
Okamoto transformations we have all 5 Dubrovin-Mazzocco solutions in [11] now derived 
via the method of -RS-pullback transformations. 

Theorem 16.21 is needed to get solutions of the equations Pyj (1/3, 1/3, 1/3, 1/5; T) and 

Pvi (1/3, 1/3, 1/2, 2/5; T) of Boalch types 37 and 43. The implied i?S'-transformations are 
T^n2( 1/3 1/2 2/5 \^^Ai^c2( 1/3 1/2 2/5 \ A. 

-"-"^4 \^ 3+3+3+1+1+1 2+2+2+2+2+2 5+5+2 y ^^^^ '"''-'4 3+3+3+I+I 2+2+2+2+2+1 5+5+1 ^ 

the end, the same normalized coverings (^i2(A^2^(a;)) and (pii{X^l{x)) can be used. But 
for intermediate computations of syzygies and application of formula (|6.12p . we may work 
with the simpler parametrized coverings ^pi2{x) and ipii{x). One convenient circumstance 
is that the point x = co does not have to be moved. 

In particular, the direct pullback of a hypergeometric equation with the local exponent 
differences 1/3, 1/2, 2/5 with respect to the covering z = ipi2{x) is a Fuchsian system with 
actual singularities at x = oo and the roots of Gi2(a;), and apparent singularities at the 
roots of Fi2, -Pi2 and Hi2. The local monodromy differences at the actual singularities 
are 4/5 or 1/3, while those differences at the apparent singularities are equal to 1 or (at 
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the roots of H12) to 2. To get rid of apparent singularities after the (imphed) Schlesinger 
transformation S{x), we have to compute syzygies between the polynomials F12, P12, -f^i2- 
To compute just a solution of Pyi (1/3, 1/3, 1/3, 1/5; T) , we apply Theorem [62] with 5 = 1, 
A = 12. The suitable syzygy is unique up to constant (in x) multiples: 

(tLi, 1, -4(3t + 4)2(2te2-t(t-8)x-(5t + 16)(2t2 + 3t + 2))), (6.18) 

where 

Li = 2t(112t^ +307t + 208)x^ - 2t(-1421t + 60t^ - 392t^ - 1040)x^ 

-2{5t + 16)(90t'* + 251t^ + 183t^ + 27t + 16)x - {2t + l){5t + 16)^(10t^ + 29t^ + 22t + 2). 

The x-root of expression (j6.12p becomes 

(2t + l)(5t + 4)(5t + 16) 
= 2(1U^ -it- 16) • ^'-^'^ 

After the normalization by conclude that a solution of Py/(l/3, 1/3, 1/3, 1/5; T37) 

is parametrized by T37 = Ai2(co) and 7/37 = Ai2(x37), like the solution in (|6.14|) . We have 
the solution yziiT^j) with T37 = r38 and 

_ 1 (2s + l)(256/ -8325^-8005*^ -3232s^ - 18445^-29505^ - 14365^-3915+64) 
^^"^ ~ 2 2w{As + 1)(6456 + 3365^ + 1045^ + 8652 - 1325 - 59) ' 

To get the parametrization obtained in 0, one has to substitute 5 ^ ~(s + l)/25. 

Similarly, the direct pullback of the same hypergeometric equation with respect to 
the covering z = ipii{x) has apparent singularities at the roots of Fn, Pn and Hn. 
For a solution of Py/ (1/3, 1/3, 1/2, 2/5; T) , we have to compute syzygies between the 
polynomials Fn, Pn, Hl^, assuming 5 = 0, A = 10. A suitable syzygy is: 

L2, x - (3t + 1)2, -30 ^*3+^^+2 ((18*2 + 13t + 9)x - (t - l)(3t + 1)^)) , (6.20) 

where L2 can be computed knowing the other two components. The x-root of expression 
()6.12p becomes 

(3t + 1)2(13*5 + 65t^ + 165*3 + 195*2 + 140t + 46) 
~ 3(3*2 + 2* + 2) (5*6 + 30*5 + 45^4 _^ 22*3 _ 13^2 _ ig* - 9) ' ^ ' ' 

After the normalization by X[i we conclude that a solution of Py7(l/3, 1/3, 1/2, 2/5; T43) 
is parametrized by T^^^ = Aii(O), 7/43 = Aii(x43). We have the solution y43(r43) with 
r43 = T42 and 

1 * + 7 

2/43 = 7: + 



2 18(*2 + 1) 

135*9 + 540*^ + 1530*^^ + 2916*^ + 3714*^ + 3486*^ + 2278*^ + 1144*2 _^ 399^ _^ 



(5*6 + 30*5 + 45^4 + 22*3 _ 13^2 _ iq^ _ 9)^3(t + 2)(* + 7) (3*2 + 3* + 2) 
To get the parametrization obtained in [5], one has to substitute * — > —(25 — l)/(5 + 2). 
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